Odd-even mass differences in self-consistent Mean-Field calculations. 
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The odd-even mass staggering (OES) of nuclei is analyzed in the context of self-consistent mean- 
field calculations. The procedure developed allows to understand the OES for spherical as well as 
for deformed nuclei. Comparison with results at the Hartree-Fock level shows the non-perturbative 
effect on this observable of the inclusion of pairing correlations. 

PACS number(s): 21.10Dr; 21.10.Hw; 21.30.-x 



Odd-even staggering of binding energies is a common 
phenomenon of several finite many-fermion systems. In 
nuclei, it has been attributed to an evidence of pairing 
correlations [Q . Assuming that masses are smooth func- 
tions of the number of neutrons and protons except for 
pairing effects, simple expressions have been derived for 
the gap parameter A based on the differences between 
binding energies of even and odd nuclei Detailed 
analyses [f| and pairing adjustments Q have been based 
on these expressions. The simplest example is the well- 
known three-point mass formula: 

(-1) N 

A {3) (N) = v ; [E(N+1) - 2E{N) + E(N-l)] (1) 

where N (Z) is the number of neutrons (protons). How- 
ever, a study of the OES in light alkali-metal clusters 
and light N=Z nuclei || has lead to the conclusion that 
this phenomenon was not due to pairing correlations but 
rather to deformation effects. That work motivated a 
study by Satula et al. JtJ about the "mean-field" contri- 
bution to the OES in nuclei, especially from deformation. 
To isolate mean-field effects, they set the pairing force to 
zero and performed Hartree-Fock (HF) calculations for 
light deformed nuclei. At this level, they already discov- 
ered an odd-even staggering, itself oscillating: 



if N is odd , 

efc-Cfc-i ^ j g even 



(2) 



where {eu — ek-i) stands for the possible gap around the 
Fermi level and is zero in spherical nuclei (apart for magic 



numbers) but different from zero for deformed nuclei be- 
cause of the spread doubly-degenerated spectrum. Thus, 
they found the deformation (Jahn- Teller effect) to be re- 
sponsible for a direct contribution to the three point odd- 
even mass formula. Consequently, in presence of pairing 
correlations A Q J d could be a measure of pairing effect 

(3) 

only, whereas A even will contain an additional contribu- 
tion responsible for the staggering of this quantity. Note 
that, such a scheme cannot account for the same oscilla- 
tion in spherical nuclei. 

The aim of the present study is to analyze for all nuclei 
the contributions to odd-even mass differences in a fully 
self-consistent mean-field picture and to extract a quan- 
tity more directly related to pairing correlations. After 
the presentation of the results, our approach will be com- 
pared to that of the Ref. . 

Finite-difference mass formulas are intended to extract 
the quickly varying part of the energy as a function of the 
neutron or proton number. The underlying assumption 
is the possible division of the energy into a smooth and 
a rapidly varying part. OES can be viewed as a probe of 
the difference of structure between odd and even nuclei, 
due to the polarization effect by the odd nucleon and to 
the reduction of pairing by its blocking. 

The energy would be smooth if those differences of 
structure would not exist. Thus the smooth part can 
be defined as ft HFBE (HFBE meaning Hartree-Fock- 
Bogolyubov Even) which is the energy obtained when all 
nuclei are calculated as if they were even ones (no block- 
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ing and no breaking of the time reversal symmetry in odd 
nuclei) . This definition of ^ HFBE has already been used 
in || in a similar context. The theoretical energy of an 
odd nucleus can be decomposed as: 



E HFB (N) = E 



HFBE 



(TV) + [E (N) — E 



pHFBE 



= E 



HFBE 



(N) + E p (N) + A(N) , 



(3) 



where E pol {N) is the difference of binding energy due 
to the polarization of the core by the odd nucleon, and 
A (A") is the positive contribution specifically related to 
the blocking effect. 

From the explicit calculation of both E HFB and 
E HFBE , one obtains EP ol (N) + A(N) through the odd- 
even mass formulas which are quantities comparable to 
experiment. The separation between E poL and A(iV) is 
not discussed in this paper. 

Following eq. H, the three point mass formula becomes: 



A 



(3) 



hfbW = A 



( hfbe(N) 



A 



(N) . (4) 

As E HFBE varies smoothly with N, one can write: 

(-i) N d 2 E HFBE 



^hfbe 



(AO 



cW 2 



(5) 



N 



We also obtain 



A 



A(iV) 



A(jV-l)+A(jV+l) 
2 



if N is odd 
if N is even, 



and similarly for A p J t (N). Similar expressions can be 
written for higher order formulas. We will use below the 
five-point formula (fourth-order): 



A 



(5) 

HFB 



(AO 



,(5) 

*HFBE 



(AT) 



,(5) 



HpoiW- (6) 



We have calculated energies using the formalism and 
forces detailed in Ref. [jlO|,[ll]]. ft is based on the self- 
consistent HFB method with an approximate particle 
number projection using the Lipkin-Nogami prescription. 
It should be mentioned that the time-odd components of 
the force are included since the time-reversal symmetry 
is explicitly broken in the calculation of odd nuclei. 

To avoid the contributions from deformation we first 
study the odd-even mass staggering along an isotopic 
chain of spherical nuclei. Seventy ground-states are cal- 
culated along the tin isotopic chain, from 100 Sn to 169 Sn. 



Each odd-N nucleus is calculated twice: first, with the 
fully self-consistent HFB scheme (several 1 quasi-particle 
(qp) configurations are investigated to get the lowest in 
energy); and second, as a HFB vacuum requiring only an 
odd average number of neutrons (HFBE case). 




A HFBE 
A 5 HFBE 



FIG. 1. Left: calculated odd-even mass differences Ag/ FB 
and AfPpg for the tin isotopic line from 100 Sn to 169 Sn. Right: 

^HFBE an0 - ^HFBE ^ 0T * ne same nuclei. 

Calculated &$ FB {N) and ^hfb(N) along this chain 
are given on the left panel of Fig. [l| A staggering is 
observed for A( 3 )(A0 whereas no staggering occurs for 
A( 5 )(A0. 

To understand these observations, let us focus first on 
the right panel of Fig. |] which displays the contributions 
of the smooth part of the energy E HFBE to A^ 3 ) and 
A' 5 ' (eq. ||and |J). Apart for the magic number A^ = 82, 
this figure shows that A^ FBE ^> A BFBE ~ 0. This 
legitimates the identification of 'E I HFBE as the smooth 
part of the energy. Indeed, the larger the order of the 
formula, the smaller the expected contribution coming 
from an existing smooth part of the energy. 

Assuming that E pol (N) + A(N) does not vary signifi- 
cantly over a few consecutive odd nuclei and gives similar 
contributions to A^ 3 ) and A^ 5 \ one can write: 

Af FB (N) - Ajj FB (N) « AP FBE (N) (7) 
AP FB (N)*A% ing+pol (N). 

To justify this equation, we compare on the left panel 
of Fig. H both sides of eq. 0. The agreement is impressive 
along the whole isotopic line (except for the magic num- 
ber A" = 82). Some one-qp configurations for odd nuclei 
are hard to converge leading to slight irregularities in the 
curve A^ 3 ) - A< 5 ) HFB. One can notice that A { BFBE {N) 
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is equal, in absolute value, for odd and even neighbors. 
The right panel of Fig. || provides a comparison between 

A hfbe and A eI p ( n ) ~ A eI p ( n ) deduced from exper- 
imentally known masses p2[ . The agreement between 
the staggering of A^ p around A^ p and &hfbe is ver y 
good as well. 
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FIG. 2. Left: A^ FBE is compared to A^ s - A 
the tin isotopic line. Right: comparison of A^ FBE 
periment (see text). 
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Eq. ^ and the above numerical results shed a new light 
on the OES problem. Indeed, the A^ 3 ) staggering must 
be understood as an oscillation around A' 5 ' due to the 
contribution to A' 3 ' of the smooth part of the energy. 
Moreover, A^ 5 ) is found to be a good measure of the 
rapidly varying part of the energy A p} iring+pol (N) which 
we want to extract while A^L still contains smooth con- 
tributions. Finally, these results fully justify the sepa- 
ration of the energy performed in eq. ||. We want to 

ir>\ 

make clear that although our calculated A HFB (N) over- 
estimates experiment (not shown here) , the oscillation of 



*HFB 



(N) is in agreement with experiment. This under- 
lines the decoupling between these two aspects. 

To extend the above analysis to deformed nuclei, we 
have calculated all cerium isotopes from 118 Ce to 166 Ce. 
The ground-state quadrupole deformation along this iso- 
topic line undergoes large variations. It goes from a re- 
gion of strong prolate deformation around 118 Ce (fa ~ 
0.37), through the spherical 140 Ce nucleus, to another 
prolate deformation region (fa w 0.31 around 160 Ce). 

On the left panel of Fig. || arc shown A H V FBE and 
A hfbe- The comparison between A^ B — A^p FB and 
A hfbe i s presented on the right panel. The same kind of 
results and agreements as for the tin isotopes is obtained: 



(3) 

HFBE 
(3) 



^ HFBE 



and A 



(3) 

HFB 



(N) - A^ FB (N) 



A { B ' FBE (N). The oscillations of A< 3 ) around A^ are 
once again well reproduced by the contribution coming 
from E HFBE . Comparison with experiment has not been 
included because the error bars on the cerium masses are 
too large to make it relevant. These calculations show 
that deformation does not modify the conclusions drawn 
above for spherical nuclei. 
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FIG. 3. Left: same as the right panel of Fig. [I] for cerium 
isotopes. Right: same as the left panel of Fig. ^| for cerium 
isotopes. 

In our separation procedure, the deformation mani- 
fests itself in several ways. First, the main effect of de- 
formation is present in Y> HFBE , apart from regions of 
shape transition and shape coexistence where the de- 
formation varies suddenly with N. The smooth contri- 
bution of the deformation to the OES is thus extracted 



through A 



(n) 

HFBE 



There is an additional contribution re- 



lated to odd-even effect. This term is usually small fyj 
and manifests itself in W° l . As seen from our numerical 
results and the above explanation, the separation given 
by eq. || is well defined whether the nuclei are deformed 
or not. Deformation will also influence the pairing part 
of the energy through correlations between particle-hole 
and particle-particle channels. The evolution of A p " iririg 
with deformation might give some informations about the 
behavior of the pairing force with respect to this degree 
of freedom. 

As it has been suggested to use HF energies as a refer- 
ence to study the A^ 3 ^ staggering , we compare in Fig. |] 
the staggering A^(2n) — A( 3 )(2n + 1) obtained with 
full HFB calculations and Hartree-Fock ones (without 
time-reversal symmetry breaking, see |?J ) for the cerium 
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isotopic chain. The single-particle splitting around the 
Fermi level (eigen-states of the HF field) (e^ — efc_i)/2 in 
even-N nuclei, due to the deformation, is also shown as it 
is presumably related to this staggering (eq. ||). The HF 
calculations do not reproduce the staggering of A^ B : 
the global trend is different between the two cases and for 
most nuclei the difference between the two calculations 
exceeds 200 keV. One can also see that the difference 
between A^ B (iV = 2n) and A < § ) FB (N = 2n + 1) is 
not a measure of the splitting of the single-particle lev- 
els (efc — efc_i)/2 around the Fermi energy. The latter 
conclusion for deformed nuclei holds of course also for 
spherical ones. 
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FIG. 4. Comparison of the staggering of A^ 3 ' for the cerium 
isotopes in HFB and Hartree-Fock calculations. The simple 
expression (efc — 6fc_i)/2 is also given. 

In summary, we have proposed a new analysis for the 
odd-even mass staggering based on self-consistent mean- 
field calculations. It assumes the definition of a "virtual" 
odd nucleus, having the structure of an even one, as the 
underlying structure of the "real" nucleus. The notion 
of underlying structure is justified by the validation of 
this energy separation obtained for a large number of 
odd nuclei, whatever the type of qp (shell, spin, parity) 
blocked self-consistently on this vacuum to describe the 
odd nuclei. Instead of being a single-particle excitation 
on an even core, an odd nucleus is seen as a qp excitation 
constructed on an odd fully paired core. It has permitted 
to identify the staggering of A^ 3 ) for spherical, as well as 
for deformed nuclei, as an oscillation around A^ 5 ) due to 
the contribution of this fully paired core. Moreover, A^ 5 ^ 
has been found to be the right quantity to extract the 
pairing part of the odd-even staggering. This conclusion 



has already been drawn from numerical comparison of 
HFBCS calculations with experiment for several possible 
theoretical quantities measuring the gap parameter || . 

Finally, our results do not confirm the ones obtained 
at the HF level that the staggering of A^ 3 ) is a measure 
of the single-particle splitting around the Fermi energy 
in deformed nuclei. A forthcoming article will present a 
full analysis of the above and complementary results. 

Methods Jl4||l5| have been recently introduced in 
which the correlations beyond mean field due to symme- 
try restorations are included. Such developments require 
a reanalysis of the different contributions to the EOS. It 
is one of the future steps of our study within the above 
framework. 
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